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Abstract 

The de Broglie-Bohm theory is about non-relativistic point-particles that move 
deterministically along trajectories. The theory reproduces the predictions of 
standard quantum theory, given that the distribution of particles over an ensem- 
ble of systems, all described by the same wavefunction ip, equals the quantum 
equilibrium distribution \ip\ 2 . Numerical simulations by Valentini and Westman 
have illustrated that non-equilibrium particle distributions may relax to quan- 
tum equilibrium after some time. Here we consider non-equilibrium distributions 
and their relaxation properties for a particular class of trajectory theories, first 
studied in detail by Deotto and Ghirardi, that are empirically equivalent to the 
de Broglie-Bohm theory in quantum equilibrium. For the examples of such theories 
that we consider, we find a speed-up of the relaxation compared to the ordinary 
de Broglie-Bohm theory. Hence non-equilibrium predictions that depend strongly 
on relaxation properties, such as those studied recently by Valentini, may vary 
for different trajectory theories. As such these theories might be experimentally 
distinguishable. 
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1 Introduction 



In the de Broglie-Bohm theory [1-3] an individual closed system is described by its 
wavefunction that satisfies the Schrodinger equation and by particle positions that move 
deterministically along trajectories with a velocity that depends on the wavefunction. 
The theory reproduces the predictions of standard quantum theory, given that the dis- 
tribution of particle positions over an ensemble of systems, all described by the same 
wavefunction if)(x, t), is given by \ip(x, t)\ 2 . The dynamics is such that if the distribution 
of particle positions equals \if)(x, to)\ 2 at a certain time to, then it equals \ip(x,t)\ 2 at 
other times t, a property called equivariance [4]. The distribution \ip\ 2 plays the role of 
an equilibrium distribution (see for example [4-6]), similar to that of thermal equilibrium 
in classical statistical mechanics, and is called the quantum equilibrium distribution. 

Current day experiments yield a very good confirmation of standard quantum theory 
and hence of quantum equilibrium. This could be accounted for by a Boltzmann-type 
argument: as shown by Diirr et al. [4] most initial configurations of the universe (relative 
to the natural measure \^(X)\ 2 dX , with \1> the wavefunction of the universe) yield the 
quantum equilibrium distribution [^(a;)] 2 for actual ensembles described by the effective 
wavefunction ip(x). In addition, it is to be expected that non-equilibrium distributions 
will tend to equilibrium and stay close to it for a reasonable time in "most" cases. 4 This 
was illustrated by Valentini and Westman [7], who presented a numerical simulation 
of such a relaxation to equilibrium for a free particle in a two-dimensional box. For a 
particular initial non-equilibrium distribution they found a rapid relaxation to quantum 
equilibrium (with relaxation time-scale roughly given by h 2 /em 1 / 2 (AE) 3 ^ 2 , where e is 
the coarse-graining length that is used in comparing the coarse-grained non-equilibrium 
and equilibrium distributions, m the mass of the particle and AE the energy spread of 
the wavefunction). 

Non-equilibrium distributions have been studied in detail by Valentini, who suggested 
possible domains where such non-equilibrium distributions might occur and how they 
might be detected, see for example [6,8-10]. In particular, Valentini suggested to look 
for non-equilibrium in astrophysical and cosmological scenarios. Assuming that the 
universe started in a special state of non-equilibrium (just like our universe seems to 
have started from a very special state of non-equilibrium on the classical statistical 
level), he considered circumstances under which this non-equilibrium might be preserved 
and even be transferred to macroscopic scales. A priori, any non-equilibrium seems 
possible. Nevertheless, in the case of a scalar field on expanding space, Valentini was 
able to come up with a quantitative prediction [11, 12]. He showed that relaxation to 
equilibrium is expected to be suppressed for a specific range of modes. One of the 

4 As explained in detail in Section 3.2, some distributions do not evolve to quantum equilibrium. For 
some wavefunctions there even is never relaxation to equilibrium. So the "most" should in the first 
place refer to some natural measure fi(dip) on Hilbert space. For example, considering product states 
il>(x\) . . . ^(xjv) (with N large), it probably holds that most non-equilibrium distributions (understood 
as empirical distributions) will tend to evolve to equilibrium, where "most" refers to the measure 
\ip(xi) . . .ip(xN)\ 2 dxi . . . dxNLi(dip). In the following we will not attempt to make these statements 
more precise. 
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possible consequences is a correction to the predictions for the temperature fluctuations 
of the cosmic microwave background in the context of inflation theory. 

In this paper we consider a certain class of trajectory theories that, just as the 
de Broglie-Bohm theory, leave \^\ 2 equivariant. In this class of theories, that was studied 
in detail by Deotto and Ghirardi [13] , the velocity field of the particles is changed by an 
additive term compared to that of the de Broglie-Bohm theory. From the point of view 
of these theories, standard quantum theory emerges as an effective theory describing 
the state of equilibrium. As such they may only be distinguished in quantum non- 
equilibrium. The goal of this paper is to examine quantum non-equilibrium and possible 
relaxation to equilibrium for these theories, by means of numerical simulations. 

In the examples of such theories that we study, the relaxation proceeds more rapidly 
than in the ordinary de Broglie-Bohm theory. This was as expected, since the additional 
term in the velocity field in general just adds to the irregularity of the motion and 
hence to its mixing-like behaviour. In particular, in the case of the ordinary de Broglie- 
Bohm theory, it was reported before that nodes in general induce chaotic behaviour 
for trajectories that come in their neighbourhood (see for example [7, 14-19]). 5 In the 
alternative theories under consideration, the additional term in the velocity fields tends 
to lead to an overall increase of its amplitude, causing the particles to travel larger 
distances, compared to the ordinary de Broglie-Bohm theory, in a given time interval. 
As such, particles tend to reach nodes sooner and undergo the chaotic motion sooner. 

Non-equilibrium predictions that depend strongly on relaxation properties, such as 
those studied recently by Valentini, may hence vary for different trajectory theories. As 
such these theories might be experimentally distinguishable. 

The outline of the paper is as follows. In Section 2, we start with reviewing this class 
of de Broglie-Bohm-type theories. In Section 3, we consider non-equilibrium distribu- 
tions and some issues that are relevant to possible relaxation. In Section 4, we present 
numerical simulations of relaxation for a number of de Broglie-Bohm-type theories, for a 
particle in a two-dimensional box (just as in the simulations of Valentini and Westman), 
finding a correlation between relaxation time and the choice of dynamics. 



2 De Broglie-Bohm theory 
2.1 Standard velocity field 

In the de Broglie-Bohm theory, an individual closed system is described by its wave- 
function ip(x,t), with x = (xi,...,xjv) € M. 3N , which satisfies the non-relativistic 
Schrodinger equation 

in^M = -J2 ^ v&(*, t) + v(xMx, t) , (i) 

k=l 



'See for example [19] or [20] for a collection of references on chaos in the de Broglie-Bohm theory. 
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and by particle positions X^i), . . . , X. N (t), whose configuration X(t) = (Xi(£), . . . ,~K N (t)) 
satisfies the guidance equation 



alt 

where the velocity field v s is given by 



<IX(f) ^v.(X(t),t), (2) 



3s (x,t) 
s,f)| J 



Vs(x,t) = i ^- 2 , (3) 



with j s = (j Sj i, . . . , j S) jv) the standard quantum probability current, given by 

j 8)fc = — Im (r V k ^) = — V fe 5|^| 2 , ^ = H exp(iS/ft) , (4) 



and which satisfies the continuity equation 



0* 



+ V-j s = (5) 



as a consequence of the Schrodinger equation. 

An arbitrary distribution p(x, t) transported along the de Broglie-Bohm trajectories 
satisfies the continuity equation 

| + v-(M = o. (6) 

Also \ip\ 2 satisfies this continuity equation, because of (5). This implies equivariance: 
if the distribution p equals \ip\ 2 at a certain time, then they are equal at all time (and 
\ip\ 2 is actually the only equivariant distribution that is also a suitably local functional 
of ip [21]). It is the generalization of the property of stationarity in classical statistical 
mechanics. 

As mentioned in the introduction \ip\ 2 plays the role of an equilibrium distribution, 
called the quantum equilibrium distribution. Given the quantum equilibrium distribution 
and the fact that measurement results are generally ultimately recorded in positions 
of macroscopic pointers, like instrument needles, computer outprint, etc., it almost 
follows immediately that the de Broglie-Bohm theory reproduces the standard quantum 
mechanical predictions (see for example [1-3,22,23]). 



2.2 Alternative velocity fields 

The flow defined by (2) and (3) is not the only one for which the distribution \ip\ 2 is 
equivariant. The continuity equation for \ip\ 2 still holds if the current j s = v s \ip\ 2 is 
replaced by 

j = 3s + ja , (7) 
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where V • j a — 0. As a result, for the velocity field 

V = W = V ' + W> (8) 

the distribution \%p\ 2 will also be equivariant. As such one can consider alternative 
theories, that just as the ordinary de Broglie-Bohm theory, agree with standard quantum 
theory in quantum equilibrium. 

This ambiguity in the choice of possible guidance equations was first studied in detail 
by Deotto and Ghirardi [13], who found that the requirement of Galilean covariance is 
insufficient to determine the guidance equation uniquely. Uniqueness could be obtained 
by extra requirements [3, 4] or by different ones [24-27]. These uniqueness results merely 
indicate the possible naturalness of certain guidance equations. 

An interesting example of this ambiguity turns up in the context of the Pauli equa- 
tion, which describes non-relativistic spin-1/2 particles [1]. Considering a single particle, 
one possible current is given by 

j s = -Im (#tvw) . (9) 

When the magnetic field is negligible and \& is a spin eigenstate, that is, \I/(x, t) = 
^(x, t)x, where x is a constant spinor, the Pauli equation will imply the non-relativistic 
Schrodinger equation for if) and the current will reduce to the one in (4). Another natural 
current is obtained by adding the divergence-free term 

ja = X (*W) (10) 

to j s . For a spin eigenstate this term survives, so that the latter current yields the 
velocity field 

+ ^W- v x W = ^ vs + 2^W V x (l*lV«x) • (u) 

This velocity field is very natural when one considers the non-relativistic spin-1/2 theory 
as the limiting case of the relativistic Dirac theory [1]. The natural current and velocity 
field for the Dirac theory namely yield the spin-term in the non-relativistic limit. The 
trajectories for this velocity field have been studied in for example [28-31]. (For other 
spins a similar ambiguity can be considered [32].) 

Note that one can have other trajectory theories in which the velocity field is not of 
the form (9), but which are still empirically equivalent with the de Broglie-Bohm theory 
in quantum equilibrium, see for example [33,34]. There even exist such theories where 
the equilibrium distribution is in general different from [■*/>[ 2 [35]. 

2.3 Some properties of the alternative velocity fields 

The velocity fields v given in (8) are not defined at the nodes of the wavefunction (that 
is, points where ip — 0). In general they also diverge near nodes. Nevertheless, in 
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the case of the ordinary de Broglie-Bohm velocity field, the set of initial configurations 
that run into a node has |-^| 2 -measure zero [36,37]. Probably similar results hold for a 
reasonable class of alternative velocity fields. 

In the case of the ordinary de Broglie-Bohm velocity field, the nodes are the only 
potential source of vorticity, while the alternative velocity fields yield vorticity even away 
from nodes. This can be seen as follows. Consider a single particle for simplicity. The 
average vorticity over a surface E with boundary <9£ is given by 

^mL^ (12) 

where § QS v • d\ is the circulation of the velocity field and A(E) is the area of the surface. 
If the surface does not contain any nodes, Stokes' theorem can be applied to write 
$8T, v ' ^1 = Js w ' d~S, where w = V x v is the vorticity vector. This vorticity vector is 
not defined at nodes (though alternatively one could define it in a distributional sense, 
see for example [38] for the case of the ordinary de Broglie-Bohm theory). 

In the case of the ordinary de Broglie-Bohm theory the velocity field is given by 
v s = V^/m, so that the corresponding vorticity vector is zero. Hence if the surface £ 
does not contain any nodes, then the average vorticity over that surface is zero. The 
average vorticity can only be non-zero when the surface contains nodes. In that case 
the vorticity is quantized [39, 40]. This is because the circulation is given by <f ds v • dl — 
ids ' dl/m = nh/m, with n integer, since S is multivalued, where different values of 
S at a point can only differ by an integer multiple of 2nh = h (since ip is single- valued) . 

Note that in the case the particle is moving in an external magnetic field B = V x A, 
the standard velocity field becomes v s = VS/m — eA/mc, so that the vorticity vector 
obtains a contribution — eB/mc. 

In the case of alternative velocity fields, the vorticity vector u> is in general different 
from zero, so that there is vorticity even away from nodes. 

In the ordinary de Broglie-Bohm theory, the chaotic motion near nodes probably finds 
its origin in the vorticity, together with the divergent behaviour of the velocity field. In 
the case of alternative velocity fields similar features hold even in the neighbourhood 
of quasi-nodes (where becomes very small). In this context, it would therefore be 
interesting to study the implications of vorticity, possibly combined with nodal or quasi- 
nodal structure, for possible chaotic behaviour. 



3 Quantum non-equilibrium and relaxation to quan- 
tum equilibrium 

3.1 Coarse-graining 

For an arbitrary distribution p one can introduce the quantity / = p/\ip\ 2 - Because p and 
| ip\ 2 satisfy the same continuity equation (with velocity field v), one has that df/dt + 
v ■ Vf = 0. This means that / is conserved along trajectories, that is, f(X(t),t) = 
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/(X(0),0) for all time t [5-8]. As such it would seem that p can never really relax to 
\ip\ 2 because their ratio is conserved along trajectories. However, relaxation should be 
understood in the sense that p goes to \ip\ 2 on a coarse-grained level [6-8]. That is, if we 
consider the coarse-grained densities p and \ip\ 2 , whose value is obtained by averaging 
over non-overlapping cells (so that their value is constant over those cells), then p/\ip\ 2 
is not necessarily conserved along the trajectories, so that p may relax to \ip\ 2 . The 
coarse-graining could be understood as corresponding to a finite accuracy of physical 
measurements. 

The deviation of p from \ip\ 2 can be quantified by means of minus the relative entropy 
of p with respect to \tp\ 2 [6-8]: 



Valentini called this quantity H in analogy with the iJ-function in classical statistical 
mechanics. It is strictly positive and zero for the quantum equilibrium distribution. 
However, since this quantity is conserved in time, just like /, it is not really suitable to 
quantify relaxation. Instead the coarse-grained if -function 



should be considered. 
3.2 Relaxation 

While we typically expect a non-equilibrium distribution to relax to equilibrium (on a 
coarse-grained level), this will definitely not happen always. That is, while we typically 
expect the coarse-grained if-function, H, to decrease over time, it may also increase 
over some time, or stay constant. This depends on both the initial distribution and the 
initial wavefunction. Let us expand on this. 

First of all, for a given wavefunction, not every initial distribution p(x) will have 
relaxed to the quantum equilibrium distribution \ip(x, T)\ 2 in a given time T. Just 
consider some distribution p(x) ^ \ip(x,T)\ 2 and evolve that distribution backwards in 
time along the de Broglie-Bohm trajectories. The resulting distribution p(x, — T) will 
then lead to the non-equilibrium distribution p(x) after time T. The distribution might 
of course relax at a later time. 

Second, for some wavefunctions there will even be no relaxation, regardless of the 
initial distribution. For example, in the case of the ordinary de Broglie-Bohm theory 
there will not be relaxation for a plane wave or a Gaussian wavefunction (see [41] for fur- 
ther examples). For a single-particle plane wave any distribution will just be translated 
over time. For a Gaussian wavefunction any distribution will either spread or contract 
over time, depending on whether the Gaussian itself is spreading and contracting. Note 
that relaxation might still occur for alternative velocity fields. 

Third, if there is relaxation of a distribution p for a particular wavefunction, then by 
time reversal invariance there exists a wavefunction for which the time evolution of the 




(13) 




(14) 
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distribution becomes time reversed [7]. The de Broglie-Bohm dynamics has the time- 
reversal symmetry (X(t),i/)(x,t)) — > (X (—t) , i/)* (x , —t)) , which reverses the direction 
in which the trajectories are traversed. As a result, an initial distribution p(x) whose 
time evolution is given by p(x,t) under the dynamics determined by ip(x,t), will evolve 
according to p(x, —t) under the dynamics determined by ip*(x, —t). 

Lastly, if the wavefunction is periodic and if the velocity field inherits this periodicity 
(which is the case for the theories for which we perform the numerical simulations), the 
trajectories typically recur. This follows by applying the standard Poincare recurrence 
theorem. One formulation of the theorem states that for a probability space (T, B, p) 
and / a measure preserving map, that is, p o f" 1 = p, one has that for A G B, for 
almost every x G A (with respect to the measure p), f n (x) G A for infinitely many 
n G N [42, p. 26]. One can apply this theorem for T C M. 3N the configuration space, 
p(dx) the measure \ip(x,0)\ 2 dx and / the map x T , where x t is the flow map associated 
to the velocity field (8) and r the period of the wavefunction. Equivariance means that 
p o x~[ 1 (dx) = \tp(x, t)\ 2 dx. Hence p o x~ l (dx) = \ip(x, r)\ 2 dx = \ip(x, 0)\ 2 dx = p(dx), so 
that the map x T is measure preserving. Note further that x t + T = x t o x T . This can be 
seen as follows. First, Y(t) = x t+T (X), X G T, is a possible trajectory, since 

^ = = v ^(m + r) = v(Y(t),t) , (15) 

with initial configuration x T (X) = X{r). Since the trajectory x t o x T (X) has the same 
initial configuration x T (X), if follows that x t + T = x t ox T . 6 Hence x™ = x nr . Application 
of the Poincare recurrence theorem now yields that for A G B, for almost every x G A, 
x nT (x) G A for infinitely many n G N, that is, the trajectory recurs an infinite number 
of times with the recurrence times given by an integer multiple of the period of the 
wavefunction. 

This recurrence theorem has the following implication for the simulations. The 
time evolution of a non-equilibrium density is simulated by considering its evolution 
only along a finite number of trajectories. From the recurrence theorem it follows that 
such a collection of trajectories will typically recur. This means that if we could run 
the simulation for arbitrary times, there would be times at which the sampled non- 
equilibrium distribution would recur. 

While the recurrence theorem indicates that distributions may evolve away from 
quantum equilibrium, it is expected that for "most" distributions the recurrence time 
is much larger than the time spent near equilibrium. 



6 Note that, except in one dimension [43], the trajectories in general do not inherit the periodicity 
of the wavefunction, that is, in general x t+nT ^ x t . Note also that the maps x t in general do not form 
a group. 
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Figure 1: The equilibrium densities at time t = 0, for the wavefunctions ip\ and ip2, 
given respectively in equations (17) and (18). 



4 Numerical simulations 

4.1 Particle in a two dimensional square box 

We consider a particle that moves in two dimensions and that is confined to a two- 
dimensional square box of side n by an infinite potential well (h, the mass m and the 
length unit are all put to 1). As such, the wavefunction is given by a superposition of 
the energy eigenfunctions 



71 



sin mil smmx2 ) 



(16) 



with energy eigenvalues E mn = (m 2 +n 2 )/2, where m and n are strictly positive integers. 

We first consider a superposition of the four lowest energy eigenstates, with equal 
weights and different phases 9 mn : 7 



2 l 

ip 1 (x 1 ,x 2 ,t) = o ( 



,i6L 



^mn(Xl,X 2 e 



m,n=l 



t 1 - 



sin mi! sin ni 2 e 



m,n=l 



(17) 

The wavefunction is periodic, with period 4n. The corresponding density \i[)i(xi,X2, 0)| 2 , 
at time t = 0, is given in figure 1. The wavefunction has only one node (within the box), 
whose trajectory is given in figures 2 and 3. 
We also consider the following wavefunction 



tp 2 (xi,X 2 ,t) 



y^e^Mxux 2 )e-^+ £ J: 



n=l v 



(j)mn(xi,X 2 )e 



(18) 



m 

m+n>2 



7 For the record, the phases are given by 9 n = 1.1525988926093297, 9 12 = 4.2775762116024665, 
6»2i = 2.1660329888555025, 6> 22 = 2.8960554218806349. They are obtained from [7], where they were 
randomly generated. 
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which also has period of An, but which doesn't have any nodes. The corresponding 
density 1^2 (an, x 2 , 0)| 2 is given in figure 1. 

In two dimensions a divergence-free two-vector field di(xi,x 2 ) can always be written 
as di = Eijdjf(xi,X2), where % is anti-symmetric with S\ 2 = 1 (as a consequence of 
Poincare's lemma), so that the velocity field of the particle can be put into the form 

/ .\ / ,\ . £ijdjf{Xi,X2,t) /in\ 

Vi{x u x 2 , t) = v 8>i {x 1 ,X2, t) + p,-— -t^- , (19) 

\lp\X\ ) X2-, tj I 

with p a constant. We will consider a number of different values for p and the following 
choices for the function /: 8 

h = M 2 , h = £ijdij s ,j , h = dij s ,i ■ (21) 

The value p = corresponds to the ordinary de Broglie-Bohm theory. The choices are 
such that there is no outgoing or incoming quantum probability flux at the boundary 
of the box (which means that the components of the currents normal to the boundary 
vanish). This means that there is no flow of \ip\ 2 in or out of the box. In figures 2 and 
3 some trajectories are given for the wavefunction ip\ given in (17). 

We consider different non-equilibrium distributions at time t — 0. The first one, 
denoted by p , is chosen to be the square of the modulus of the ground state, that is, 

p (x 1 ,x 2 ,0) = \(pu(x 1 ,x 2 )\ 2 = sin 2 x x sin 2 x 2 ■ (22) 

The other non-equilibrium distributions, denoted by pi, % — 1, . . . , 4, are obtained from 
p by translation and contraction, that is, 

p 1 (x 1 ,x 2 , 0) = 4A)(2a;i, 2x 2 , 0) , (x l7 x 2 ) e [0, tt/2] x [0, tt/2] , 

P2(x 1: x 2 , 0) = 4p (2(x 1 - vr/2), 2x 2 , 0) , (an, x 2 ) G [tt/2, tt] x [0, tt/2] , 

p 3 ^i, x 2 , 0) = 4p (2x!, 2(x 2 - tt/2), 0) , (an, x 2 ) e [0, tt/2] x [tt/2, tt] , 

p 4 (xi, x 2 , 0) = 4p (2(an - tt/2), 2(a; 2 - tt/2), 0) , (an, x 2 ) e [tt/2, tt] x [tt/2, tt] , (23) 

and are zero outside the domains specified. So these distributions have their support in 
the four different quadrants of the square box. 



8 Notc that the velocity field (11), for the case of a spin-eigenstate \ = (1 0) T an d for a factorizable 
wavefunction ip(xi, x 2 , xz,t) — (p(xi,X2,t)t](x3,t), reduces to 

vi = -d 1 s v + — ^,d 2 M 2 , v 2 = -d 2 s v -- h r ^dM\ v 3 = -d 3 s TI , (20) 

where S v and S v are respectively the phases of ip and 77. Hence in this case the motion in the (xi, x 2 )- 
planc decouples from that in the X3 direction and the velocity field in the (xi, X2)-plane is of the form 
(19), with pt = h/2m and / = /1 = \ip\ 2 . 
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Figure 2: Figure (a) displays the trajectory of the node of the wavefunction ipi. The 
other figures display particle trajectories for different values of /i and /, for the wave- 
function ip\. In each case, the initial position is the center of the box. The final position 
is denoted by an arrow. All trajectories are for the time interval [0, 4tt]. Note that the 
trajectory in figure (e) does not display periodicity over this time interval. 
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Figure 3: Space-time trajectories of the node (bold lines) and particle trajectories (thin 
lines), for the wavefunction ipi, over the time interval [0, 47r]. In each case, the initial 
position is the center of the box. 
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4.2 Some details of the algorithm 



The algorithm to calculate the evolution of the non-equilibrium densities is similar to 
that of Valentini and Westman [7]. Just as in their work, the calculation of the trajec- 
tories is based on the Runge-Kutta-Fehlberg algorithm with Cash-Karp parameters [44] 
(but with a starting value of A equal to 1CT 5 and the maximal number of time steps 
equal to 10 5 ). 

In order to calculate the time evolved non-equilibrium density p(x\,X2,t) from the 
initial density p(x±,X2,0), a uniform lattice is placed over the box, whose lattice points 
have the coordinates (A;tt/1024 - tt/2048, Ztt/1024 - vr/2048), where k, I = 1, 2, . . . , 1024 
(with hence a total of 1024 x 1024 lattice points). Each lattice point is then evolved 
backwards in time from time t to time t — using the de Broglie-Bohm dynamics. 
Then, using the constancy of / = p/\ip\ 2 along a trajectory (see Section 3.1), one has 
that the time evolved density at a lattice point with coordinates (xi,X2) is given by 
p(xi,x 2 ,t) = \ip(xi,X2,t)\ 2 f(xi(0), #2(0), 0), where (xi(0), #2(0)) are the coordinates of 
the backtracked position. This method of using the backtracked positions is advanta- 
geous compared to one where positions are evolved forward in time [7]. If the calculation 
of the backtracked position involves more than the 10 5 time steps, it is halted and the 
lattice point is ignored when calculating the coarse-grained and smoothed density. 

The coarse-graining of a density p is done by averaging over square non-overlapping 
cells of side 7r/32 (with a total of 32 x 32 cells). Each such coarse-graining cell C 
contains 1024 lattice points, of which a certain number N c can be backtracked. So the 
coarse-grained density p in a certain cell C is given by J2 C p(xi,x 2 ,t)/N c , where the 
sum ranges over coordinates (x±,X2) of lattice points that are contained in the cell C 
and that could be backtracked. 

Lattice points near nodes are generally harder to backtrack because the velocity 
field generally diverges there. This is in particular the case for the lattice points near 
the boundary of the box. Therefore in order to speed-up the calculation of the time 
evolved non-equilibrium densities, we have ignored those lattice points that lie within 
two coarse-graining cells near the boundary. As such there is an error in calculating 
the coarse-grained if-function, since the ignored regions near the boundary yield no 
contribution to the integral. Therefore our calculated coarse-grained if-function only 
measures the difference between the calculated distributions p and \ip\ 2 over the region 
within the box that excludes the ignored coarse-graining cells. 

In the figures, smoothed densities p are displayed, instead of the coarse-grained 
densities p. The reason for the smoothing is to average out fine-grained fluctuations. 
It is done by averaging over square cells of side 11/ 16, ignoring again the cells near 
the boundary. The cells may be overlapping and can be mapped to one another by 
translations over an integer multiple of 7r/108 along the x\ and x 2 directions. More 
explicitly, we consider the points with coordinates (&7r/128 + 37r/32, Itt/128 + 37r/32), 
with k, I = 0, 1, 2, ... , 104, and calculate p at those points by averaging the density p 
over cells of side 7r/16, ignoring again the lattice points that could not be backtracked 
(so that each average is done over at most 4096 lattice points). In generating the plots 
with Matlab an interpolation of those densities values is then introduced. 
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* = 8tt 






t = 2tt 


* = 4tt 




t = 4tt 


p = 





99.95 


99.89 


n = 





99.89 


99.62 


p = 


h f 2 


100. 


99.62 


p = 


i, h 


99.98 


99.52 


p = 


1, / 2 


99.99 


91.71 


p = 


2,/i 


99.52 


94.44 


n = 


1, /s 


100. 


99.82 










p = 


2, /a 


99.99 


97.16 



Table 1: Percentage of backtracked lattice points, for the wavefunction ipi, and for the 
guidance equations with different values of p and / (up to two decimal places). 





t = Qn t=12n 


p 


= 


100. 100. 


p 


= 2, h 


99.45 89.21 



Table 2: Percentage of backtracked lattice points, for the nodeless wavefunction ip 2 , for 
the ordinary de Broglie-Bohm guidance equation and for the one with p — 2 and / = f\ 
(up to two decimal places). 

4.3 Results of the numerical simulations 

Figures 4-7 display the time evolution of the distributions l^il 2 an d of the different non- 
equilibrium distributions p^, i — 0, . . . , 4 for ipi (given by (17)), and this for different 
guidance equations (19) that are distinguished by different values of p and /. Figure 8 
displays the distributions for the nodeless wavefunction ip 2 (given in (18)). 

The percentages of lattice points that could be backtracked in each case are given 
in the tables 1 and 2. In most cases the percentages over coarse-graining cells did not 
differ all too much from the percentages over the box. The coarse-graining cells with 
the worst percentages were generally located near the boundary of the box. 

The calculated values of the coarse-grained ^/-function are given in tables 3 and 
4. Note that some values of H are negative, which is theoretically impossible. In 
general there will be some deviation from the actual values, because of several reasons. 
First, these densities are calculated using only a lattice sampling, second there are the 
numerical errors that come up in calculating the time evolution of the distributions, and 
third, by ignoring the contributions of the regions near the boundary in the integral that 
defines H, we have effectively put p zero in those regions. In [7] Valentini and Westman 
estimated the error in H to be about 2%. Assuming a similar error in our case, we have 
presented the values of H rounded to the nearest integer. 

From the pictures, as well as from the values of the coarse-grained if-function, it 
is clear that the non-equilibrium distributions evolve closer to equilibrium over time. 
Furthermore, the relaxation time decreases when the value of p increases. As mentioned 
before, this might be explained by the fact that the extra term in the velocity adds to 
the overall irregularity of the motion and to the vorticity in particular. 
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t=0 t=4re t=8jt: 




x 2 x i x 2 x i x 2 x i 



Figure 4: Smoothed densities |-?/>i| 2 and p at times t = 0, t = Air and t = 8tt, for guidance 
equations with different values of \x and / = f\, for the wavefunction ipx. Notice that 
the relaxation to equilibrium is better for higher value of fi. 
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t=0 



t=4jt 



t=8it 




Figure 5: Smoothed densities p% and p 2 at times t = 0, t = Air and t = 8ir, for guidance 
equations with different values of fi and f — fi, for the wavefunction ipi. 
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Figure 6: Smoothed densities ps and p4 at times t = 0, t = An and t = 8ir, for guidance 
equations with different values of p and f = fi, for the wavefunction ipi. (In the case 
of p4, p — and t = 8tt, one peak in the distribution is not fully plotted.) 
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Figure 7: Smoothed densities \ipi\ 2 and po at times t = 0, t = 2tt and t = Air, for 
guidance equations with different values of /i and /, for the wavefunction (In the 
case of po, \i = and t = 2ir, one peak in the distribution is not fully plotted.) 
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t=0 t=6:t t=12st 




Figure 8: Smoothed densities \ip2\ 2 and po at times t = 0, t = 6tt and t = 12n, for the 
ordinary de Broglie-Bohm guidance equation and for the one with \i = 2 and f = fi, for 
the wavefunction 
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54 
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Po, P = 
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23 


21 
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fi 


169 


26 
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Po, P = 


l,/2 


54 


22 


18 


P2, P = 


2, 


fi 


169 


17 


1 


Po, P = 


l,/3 


54 


5 


1 


P~3, P = 







130 


94 


61 


Po, P = 


2,/ 3 


54 


2 


-1 


P~3, P = 


1, 


fi 


130 


12 


1 












P~3, P = 


2, 


fi 


130 


5 















Pi, P = 







300 


207 


158 












Pi, P = 


1, 


fi 


300 


54 


19 












Pi, P = 


2. 


fi 


300 


21 


7 













Table 3: Values of H (rounded to the nearest integer) for the densities plotted in figures 
4-7. Notice the decrease of H with increase of /i. 





t = 


t = 6vr 


t = 12tt 


Po P = 





17 


17 


16 


Po P = 


2, fi 


17 


1 


-1 



Table 4: Values of H (rounded to the nearest integer) for the densities plotted in figure 
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Note that in the case of the ordinary de Broglie-Bohm velocity field, the relaxation is 
not as good as in the simulation of Valentini and Westman [7] , where a superposition of 
the first 16 modes was considered instead of the first 4. This is probably related to the 
fact that more chaotic behaviour is expected when the number of nodes increases [14]. 
The wavefunction of Valentini and Westman namely has much more than one node. 

In the case of the guidance equation with / = f 2 there is not a very good relaxation 
neither for \i = 0.5, nor for fx — 1, compared to the other alternative guidance equation 
equations (though it is still better than for the ordinary de Broglie-Bohm theory). We 
have tried to perform the simulations also for \x = 1.5, but found that only 76.78% of 
the lattice points could be backtracked. Since this percentage is rather low, we did not 
include these results in the paper. Figures 2 and 3 display some trajectories for different 
possible guidance equations. The trajectory for / = f 2 and /i — 1 circles around the 
node from time t — to t — 47r, unlike the trajectories corresponding to the other 
guidance equations, which tend to cover a bigger area of the box in that time. This 
might be a possible reason why the relaxation is not so good in this case. If the dynamics 
is such that initial positions that start near a node remain close to it for long periods of 
time and merely circle around it, one might expect that there is insufficient mixing-like 
behaviour for relaxation to occur efficiently 

In the case of the wavefunction without node, there is some convergence towards 
quantum equilibrium for the ordinary de Broglie-Bohm guidance equation. However, 
there is definitely not a good relaxation. It is unclear how much better the convergence 
would be for larger times. In the case of the guidance equation with /i — 2 and / = 
fi there is an approximate relaxation to equilibrium. Presumably the more efficient 
relaxation in the latter case is related to possible chaotic behaviour induced by the 
vorticity. In the ordinary de Broglie-Bohm theory there is probably no chaotic motion 
for this wavefunction, due to the absence of nodes. 

5 Conclusion 

We have considered non-equilibrium distributions and their relaxation properties for 
a class of de Broglie-Bohm-type theories, and found that the relaxation time depends 
substantially on the form of the guidance equation. In particular, in the case of the 
examples of such theories considered here, there was a faster relaxation compared to the 
ordinary de Broglie-Bohm theory. 

Only non-relativistic quantum theory was considered, but similar results are expected 
to hold for relativistic quantum theory or quantum field theory. As such, these results 
might have potential implications for empirical predictions concerning non-equilibrium 
distributions, like for example those of Valentini [11, 12] that were mentioned in the 
introduction. But, with this being said, Valentini's predictions involve a field ontology 
for a scalar field, with a velocity field that is most natural. In support of his choice, 
Valentini actually proves uniqueness of this velocity field for the vacuum state (the state 
of interest for his analysis) under certain natural assumptions [12]. 

Actually, when it comes to quantum field theory, there is, apart from the ambiguities 
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in the velocity field, also a potential ambiguity in the choice of ontology. Some de Broglie- 
Bohm-type models for example suggest a field ontology, while others suggest a particle 
ontology (see for example [45,46] for possible particle ontologies and [47] for a review on 
possible field ontologies). Such ambiguities may potentially be relevant too for possible 
empirical predictions concerning non-equilibrium. 

Note that similar ambiguities arise in for example the context of Nelson's stochas- 
tic mechanics [48,49]. Nelson's mechanics [50,51] can be regarded as a modification 
of the de Broglie-Bohm theory in which the deterministic de Broglie-Bohm dynamics is 
replaced by a diffusion process (though Nelson himself wanted to derive the wavefunc- 
tion and Schrodinger's dynamics from the stochastic process instead of assuming them 
independently). This diffusion process is expected to yield convergence to equilibrium 
for most distributions (see [52] for some rigorous results). As with the de Broglie-Bohm 
theory there are many ways in which a Nelsonian type of dynamics can be introduced 
that preserves the equilibrium distribution \ip\ 2 . In particular, there are Nelsonian the- 
ories with arbitrary values of the diffusion constant [48]. In this case, relaxation time is 
expected to decrease with increasing diffusion constant. 
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